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ABSTRACT 


Transplutonium isotopes are scarce and need to be produced by irradiation in 
high flux reactors, but their production is inefficient and optimization studies are 


needed. This paper analyzes the physics nature of transplutonium isotopes production 


28 
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by taking Cf-252, Cm-244, Cm-242 and Pu-238 as examples. Traditional methods 
based on the Monte Carlo burnup calculation are faced with the shortcomings of large 
amounts of calculation and are unable to analyze the individual energy intervals in 
more detail, thus cannot support the refined evaluation, screening and optimization of 
the irradiation schemes. After grasping the physics nature and simplifying the 
complexity of the production process, we proposed a rapid diagnosis method for 
evaluating the radiation schemes based on the concept of “Single Energy Interval 
Value (SEIV)" and “Energy Spectrum Total Value (ESTV)". The rapid diagnosis 
method can not only avoid the tedious burnup calculation but also help to provide 
direction for optimization. The optimal irradiation schemes for producing Cf-252, 
Cm-244, Cm-242 and Pu-238 are determined based on the rapid diagnosis method. 
The optimal irradiation schemes can greatly improve production efficiency. Compared 
with the initial scheme, the optimal scheme improves the production efficiency of Pu- 
238 by 7.41 times, 11.98 times for Cm-242, 65.20 times for Cm-244 and 15.08 times 
for Cf-252, respectively. The work in this paper realizes the refined analysis of 
transplutonium isotopes production and provides a theoretical basis for improving 


production efficiency. 


Keywords:  lransplutonium isotopes, Rapid Diagnosis Method, Production 


optimization, Single Energy Interval Value, Energy Spectrum Total Value. 


1. INTRODUCTION 
Transplutonium isotopes [1] refer to nuclides with an atomic number greater than 
or equal to 94 (plutonium element). Most are radioactive isotopes with unstable nuclei 
and can spontaneously release rays (a, f, y) or neutrons by spontaneous fission. 


Transplutonium isotopes mainly include: the 94" element plutonium (Pu), the 95" 
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element americium (Am), the 96" element 96 curium (Cm), the 97" element 
berkelium (Bk), the 98^ element californium (Cf), et al. All transplutonium isotopes 
are metals, as shown in Figure 1. Because of their research and application value, 


transplutonium isotopes have attracted extensive scientific attention [2-4]. 


Figure 1. The photographs of plutonium, americium and berkelium, californium 

Transplutonium isotopes are more scarce and expensive than medical isotopes 
[5-8]. For example, Cf-252 is the most expensive element, at $27m a gram, 650,000 
times more than gold. Cf-252 is a strong neutron source, emitting 2.35x10* neutrons 
per microgram per second with a half-life of 2.65 years. Cf-252 releases neutrons in a 
spectrum similar to a fission reactor, making it an ideal source of neutrons for reactor 
start-up [9]. Cf-252 is also effective in the treatment of cervical cancer [10], in oil 
exploration [11], and in prompt gamma neutron activation analysis (PGNAA) [12], et 
al. Pu-238 is widely used in the radioisotope thermoelectric generator because of its 
moderate half-life, high thermal power density and easy radiation protection from a 
decay [13]. Cm-242 and Cm-244 are portable a-ray sources that can power satellites 
and spacecraft, and they are also used as target materials to produce Cf-252 [14]. 

Most transplutonium isotopes do not exist naturally and must be produced [15- 
18]. There are three methods of production: 1) production in accelerators [19], 2) 
production. during a thermonuclear explosion [20], and 3) production during 
irradiation in a high flux reactor [21]. Among them, irradiation in a high flux reactor 
is the most stable and efficient production method, and it is the only way to produce 
Cf-252 in commercial operation. Only two laboratories worldwide can stably produce 
Cf-252 using a high-flux irradiated curium target [22]. One is the Oak Ridge National 


Laboratory (ORNL) in the United States, and the other is the Research Institute of 
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Atomic Reactors (RIAR) in Russia. However, worldwide production of 
transplutonium isotopes remains low. In 78 production runs over the past 60 years 
[23], only 1.2 grams of Bk-249, 10.2 grams of Cf-252, 39 milligrams of Es-253, and 
15 picograms of Fm-257 were produced through the High Flux Isotope Reactor 
(HFIR) [24] and the Radiochemical Engineering Development Center (REDC) [25] in 
the United States. 

Neutron capture reactions and decay reactions occur during the production by 
reactor irradiation [26-27]. For each neutron absorbed by a nuclide, the atomic 
number of the nuclide increases by one, becoming a new nuclide. So the target 
material changes from nuclides with a small mass number to nuclides with a large 
mass number through a series of nuclear reactions (radiation capture, f-decay). Figure 


2 shows the nuclide conversion process of producing transplutonium isotopes in a 


high flux reactor [28]. 
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Figure 2. Nuclide conversion process of producing transplutonium isotopes 


Besides the absorption reaction, the fission reaction also occurs. Once the fission 
reaction occurs, the production of transplutonium isotopes is stopped, causing the 


fission losses and leading to a low conversion rate. Therefore, the fission reaction 


should be reduced while the absorption reaction should be increased to improve 
production efficiency. Because the cross section is related to the neutron spectrum 
[29-30], the neutron spectrum around the target should be studied. The neutron 
spectrum around the target greatly influences production efficiency [31]. To 
summarize the literature review, two questions should be answered to improve the 
production efficiency of transplutonium isotopes: 1) how to find the optimal neutron 
spectrum, and 2) how to achieve the optimal neutron spectrum. The background of 


optimization of transplutonium isotopes production is shown in Figure 3. 
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Figure 3. The background of transplutonium isotopes production optimization 


For answering the first question, it is difficult to find the optimal neutron 
spectrum because of the high complexity of the production process, and there is still 
no corresponding Neutronics model that can analyze each energy interval 
individually. Traditional methods search for the optimal neutron spectrum by a large 
number of Monte Carlo burnup calculations. The burnup equation is independent of 
energy, so it is impossible to perform multi-group burnup calculation. Therefore, the 
traditional methods take many computations, and can only macroscopically analyze 
the whole energy spectrum, not the particular energy intervals in more detail, thus 


cannot support the refined evaluation, screening and optimization of the irradiation 
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scheme. So we need a rapid diagnosis method for transplutonium production. 
Grasping the essence of the production process and simplifying the complexity to 
achieve the refine Neutronics analysis, so as to support the optimization, improve the 
production efficiency and reduce the production cost. 

For answering the second question, the neutron spectrum has a complex 
relationship with the design parameters, so it is necessary to perform reactor physics 
calculations [32-34] and sensitivity analysis. After the sensitivity analysis of some 
design parameters to production efficiency, we can achieve the optimal neutron 
spectrum by optimizing these design parameters, such as 1) the selection of the 
moderator, 2) the layout of the moderator, 3) the position of the target, 4) the size of 
the target, and 5) the shape of the target. We try to adjust these parameters to obtain a 
production scheme with high efficiency, that is, by adjusting these parameters to 
obtain an optimal neutron spectrum for producing transplutonium isotopes. 

This paper analyzes the optimization method of transplutonium isotopes 
production by taking Cf-252, Cm-244, Cm-242 and Pu-238 as examples. By 
proposing the concepts of “Single Energy Interval Value (SEIV)” and “Energy 
Spectrum Total Value (ESTV)", we establish a rapid diagnosis method and realize the 
detailed analysis of the energy spectrum around the target. The rapid diagnosis 
method can not only avoid the tedious burnup calculation but also help to provide 
direction for optimizing the irradiation scheme. Then, an optimal irradiation scheme is 
proposed after the sensitivity analysis based on the rapid diagnosis method. The paper 
is organized as follows: Chapter 2 introduces the rapid diagnosis method, Chapter 3 
gives numerical and experimental verification, Chapter 4 introduces the optimal 


process based on sensitivity analysis, and Chapter 5 gives a conclusion. 
2. RAPID DIAGNOSIS METHOD 


2.1 The High Flux Reactor 


All analyses are performed in a high flux lead-bismuth reactor that is under 


design. It has a 90-day refueling cycle with the highest flux of 6.7x10'°. The inlet 
temperature is 170°C the outlet temperature is 536.5°C and the coolant velocity is 4.0 


m/s. ?*Pb-Bi is used as the coolant, and "*Pb is used as the reflector layer. Figure 4 
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4 shows the X-Y and X-Z sections of the reactor. 
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L3 


5 VIENT 
6 Figure 4. The X-Y and X-Z sections of the reactor 
T The fuel rod is a simplified model consisting of fuel pellets, air gap and cladding, 


8 with a total length of 50 cm. The fuel assembly has a hexagonal design and contains 


9 91 fuel rods. The geometry of the fuel rod and the fuel assembly are shown in Figure 


10 5. The detailed parameters of the high flux reactor are given in Table 1. 


Fuel rod Fuel assembly 
11 
12 Figure 5. The geometry of the fuel rod and the fuel assembly 
18 Table 1. The detailed parameters of the high flux reactor 
Parameters Value 
Thermal power 150 MW 
Loading capacity of Fuel / U-235 779 [175.3 kg 
Equivalent diameter of active zone 58.14 cm 
Height of active zone 50 cm 


Average line power density 302.45 W/cm 


Maximum line power density 404.20 W/cm 
Average volume power density 1130 W/cm? 
Fuel material U-10%Zr alloy 
Enrichment of ^*U 25 96 

Inner / outer diameter of fuel rod 4/ 4.6 mm 
Filled gas in the gaps Helium 

Fuel rod clearance width 0.1 mm 
Cladding material T91 
Thickness of the cladding 0.2 mm 
Ratio of pitch to diameter (P/D) 1.1304 
Number of assembly 109 

Number of fuel rods per assembly 91 

Axial / radial thickness of reflector layers 80 / 120 cm 


The target for producing Cf-252 is a mixture of plutonium, americium and 
curium [9], and the nuclide composition of the target is shown in Table 2. The target is 
irradiated for five cycles, and each cycle consists of a 25-day irradiation. This target 
irradiation is representative of typical operations during a californium production 
campaign. 


Table 2. The nuclide composition of the target for producing Cf-252 


Nuclide Weight Nuclide Weight 
Pu-238 0.002 g Cm-244 8.734g 
Pu-239 0.001 g Cm-245 0.271g 
Pu-240 0.565 g Cm-246 22.349g 
Pu-242 0.016 g Cm-247 0.647g 
Am-243 3.903g Cm-248 4.743g 


2.2 Production Evaluation by Monte Carlo Burnup Calculation 

The yields of transplutonium isotopes can be obtained by performing the Monte 
Carlo burnup calculation, which is the coupling of the Monte Carlo criticality 
calculation and the burnup calculation. Figure 6 gives the flowchart for evaluating the 


transplutonium isotopes production by the Monte Carlo burnup calculation. 
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Figure 6. Flowchart for evaluating production by Monte Carlo burnup calculation. 
The point burnup equation describes the nuclides’ transmutation with time when 
the target is irradiated in a high-flux reactor. For each nuclide in the burnup chain, its 


time-dependent point burnup equation can be written as 


; (1) 


where n; is the density of the i" nuclide, A47" is the effective decay constant of the ;^ 


nuclide, and 5;/" is the branching ratio for transmutation of the 7" nuclide to the ^ 


nuclide. 4" and b, can be calculated from the following formula, 


; (2 
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where 4; is the decay constant of the ;^ nuclide, ¢ is the neutron flux, and o;; is the 
y J 


one-group cross section where the i" nuclide's reaction generates the j^ nuclide. 
During the burnup calculation, the one-group cross section and neutron flux 

within each burnup step are assumed to be constant so that the point burnup equation 

can be treated as a first-order linear ordinary differential equation. Writing it into a 


matrix form, we have the following more concise expression, 


; (3) 


where matrix A represents the coefficient matrix of the N order point burnup equation, 
which is regarded as a constant matrix in a single burnup step. 
Given the initial boundary conditions, the solution of Equation (3) can be written 


in the form of an exponential matrix, 


; (4) 


where matrix At is referred to as the burnup matrix for the convenience of quotation. 

Because of the extremely low yields of transplutonium isotopes and the complex 
nuclide composition of the target, the point burnup model usually adopts a complex 
nuclide system, which needs to cover more than a thousand isotopes, including a large 
number of short half-life nuclides. For example, the ORIGEN-S burnup database [35] 
contains about 1,500 nuclides. Some of these nuclides, such as Po-212, have a half- 
life of 107s. Therefore, the point burnup equations of transplutonium isotopes 
production are not only large in number but also rigid, which brings difficulties to its 
solution. 

As shown in Equation (2), the cross section of the target is required for the 
Monte Carlo burnup equation, so the criticality calculation should be performed. The 


k-eigenvalue neutron transport equation can be written as [36] 


; (5) 


where L represents the leakage operator, C represents the collision operator, S 


27 


28 


represents the scattering operator, M represents the fission operator, keff represents the 
effective multiplication factor, and ó is the neutron flux. 

Physical parameters around the target, such as the neutron flux, fission reaction 
rate and absorption reaction rate, can be obtained by solving Equation (5). The Monte 


Carlo method has the following formula for calculating reaction rates, 


(6) 
where Y, represents the macroscopic cross section in the g^ energy interval, and > 
represents nuclides’ cross section related to position r and energy E. 

Therefore, we have the following formula for calculating the cross section 

needed by Monte Carlo burnup calculation, 
(7) 

The above formulas show that the Monte Carlo burnup calculation is complex 
and costly [37-38]. Tested in the high flux reactor shown in Figure 4, with the 
calculation parameters of 1,000,000 neutrons per cycle, 100 inactive cycles, 300 
active cycles and 10 burnup steps, it takes 150 minutes to perform one Monte Carlo 
burnup calculation with 64 threads on the 2nd Gen AMD EPYC™ 7742. Even though 
only considering the production of Cf-252, Cm-244, Cm-242 and Pu-238, and only 
analyzing the five design parameters (shown in Figure 9), a total of 9720 calculations 
are needed, with a computing time of 1458,000 minutes (1012.5 days). This amount 
of computation time is unacceptable. 

Moreover, since the Monte Carlo burnup calculation only needs the one-group 
cross section, the calculations can only give the influence of the whole neutron 
spectrum on the production efficiency, not a closer analysis of each single energy 
interval. So when optimizing the neutron spectrum, we can only from the macroscopic 
design parameters, cannot know which energy intervals are favorable and which are 
unfavorable to production, and thus cannot provide more detailed guidance for the 
optimization. 


Therefore, we need to grasp the physical nature of the production process and 


find a new rapid diagnosis method which can simplify the calculation and analyze the 
production process in more detail. A rapid diagnosis method is helpful for evaluating, 


screening and optimizing the irradiation scheme. 


2.3 Rapid Diagnosis Method Based on SEIV and ESTV 

The root cause of the low production of transplutonium isotopes is fission loss 
[14]. When the target is irradiated in a high flux reactor, the target material changes 
from nuclides with a small mass number to nuclides with a large mass number 
through absorption reactions. Besides the absorption reaction, the fission reaction also 
occurs. Once the fission reaction occurs, the production of transplutonium isotopes is 
stopped, causing fission losses and leading to a low conversion rate. Therefore, the 
ratio of absorption rate to fission rate (A/F) is an important physical quantity affecting 
production efficiency. The production efficiency can be assessed by analyzing the 
A/Fs of the nuclides in the nuclide chain. Moreover, the A/F can be calculated in 
individual energy intervals, thus achieving refined Neutronics analysis. Taking the 
production of Cf-252 from the target described in Table 2 as an example, Figure 7 
shows the A/Fs of the intermediate nuclides, Cm-244, Cm-245, Cm-246, Cm-247, 


Cm-248, Bk-249, Cf-250 and Cf-251. The nuclear data comes from JANIS [39]. 
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Figure 7. The A/Fs of the Intermediate nuclides for producing Cf-252 

As shown in Figure 7, the A/Fs of different nuclides vary greatly, and so do the 
A/Fs of different energy intervals, so different nuclides and different energy intervals 
have significantly different effects on production efficiency. It is not accurate to only 
analyze the whole energy spectrum, and a more detailed analysis of each energy 
interval is needed. 

A rapid diagnosis method based on the neutron spectrum around the target is 
proposed to realize the rapid and refined evaluation of the different production 
schemes. The model only needs to obtain the neutron spectrum around the target by 
the Monte Carlo criticality calculation [40-41] to evaluate the efficiency of different 
production schemes, avoiding the burnup calculation and greatly reducing the 
required computing resources. Therefore, the rapid diagnosis method can screen out 
the optimal scheme for transplutonium isotopes production quickly, and helps provide 
direction for optimizing the irradiation scheme. 


The model is divided into four steps: 1) the neutron spectrum around the target is 


28 


29 


obtained after the Monte Carlo criticality calculation, which can be described by the 
46-group neutron flux, the 46-group fission rate and the 46-group absorption rate, 2) 
define the physical quantity “Single Energy Interval Value (SEIV)", and calculate the 
SEIV of these 46 energy intervals, 3) define the physical quantity “Energy Spectrum 
Total Value (ESTV)", and use the 46-group neutron flux and the 46-group SEIV to 
calculate the ESTV of the neutron spectrum, and 4) evaluate the efficiency of different 
schemes by comparing the ESTV. Those with higher ESTV are considered to have 
higher production efficiency and yield. The schematic diagram is shown in Figure 8. 


The Single Energy Interval Value (SEIV) is calculated by 
; (8) 


represents the serial number of energy intervals, Ra; represents 


ee 5*5 
l 


where the subscript 
the absorption rate in the i” energy interval, Aj; represents the fission rate in the ;^ 
energy interval, ó; represents the neutron flux in the 7" energy interval. 

The SEIV is the product of two expressions. The first is the A/F, indicating 
whether it tends to undergo an absorption or a fission reaction within this energy 
interval. The other is the ratio of the absorption rate to the neutron flux, which 
indicates the probability of the absorption reaction occurring in this energy interval. 
Therefore, a higher SEIV means that within this energy interval, absorption reactions 
tend to occur rather than fission reactions, and more absorption reactions occur. The 
energy interval with a high SEIV is considered important in the production process, 


and the neutron flux in this energy interval should be increased as much as possible. 


The Energy Spectrum Total Value (ESTV) is calculated by 
(9) 


The ESTV is the integral value of SEIV along with the neutron spectrum, 
indicating the target nuclide's overall trend of absorbing neutrons without fission. We 
believe that the production scheme with higher ESTV produces transplutonium 
isotopes more efficiently. The ESTV 1s only related to the neutron spectrum around 


the target, and the burnup calculation is no more needed, so the ESTV can be used for 
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the rapid evaluation of the production scheme. 
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Figure 8. The schematic diagram of the rapid diagnosis method 


3. NUMERICAL AND EXPERIMENTAL VERIFICATION 


3.1 Numerical Verification 

The rapid diagnosis method is verified by comparing the burnup calculation with 
the ESTV. We optimize the production efficiency by adjusting the following five 
design parameters: 1) the selection of the moderator, 2) the layout of the moderator, 3) 
the position of the target, 4) the size of the target, and 5) the shape of the target. The 


optional values for these five design parameters are shown in Figure 9. 
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Figure 9. The optional values for these five design parameters 

We use a large hexagonal prism target with Position #1 as the initial scheme, and 
zirconium hydride is selected as the moderator. Different irradiation schemes are 
obtained by modifying the layout of the moderator. We perform criticality calculations 
to obtain the neutron spectrum around the target. The SEIV and ESTV are determined 
and used to evaluate the efficiency of the irradiation schemes. 

We take Pu-238, Cm-242, Cm-244, and Cf-252 as examples. The target for 
producing Pu-238 is Np-237. The target for producing Cm-242 is Am-241. The target 
for producing Cm-244 is Pu-239. The nuclide composition of the target for producing 
Cf-252 is shown in Table 2. The SEIV and ESTV for the target without the moderator 
(bare target) and the six layouts of the moderator are calculated, so there are seven 
production schemes for each transplutonium isotope to be produced. The SEIV for 


producing these four transplutonium isotopes is given in Table 3. The ESTV for 


producing these four transplutonium isotopes is given in Table 4. 

As shown in Table 3, the optimal neutron spectrum for producing Pu-238, Cm- 
242, Cm-244 and Cf-252 is not always the softer, the better. The energy interval 
[4.14E-07 MeV, 8.76E-06 MeV] has the highest SEIV for producing Pu-238. The 
energy interval [1E-11 MeV, 1E-07 MeV] has the highest SEIV for producing Cm- 
242. The energy interval [1E-07 MeV, 4.14E-07 MeV] has the highest SEIV for 
producing Cm-244. The energy interval [8.76E-07, 1.86E-06] has the highest SEIV 
for producing Cf-252. Moreover, the SEIV varies greatly in different energy intervals. 
For example, the maximum deviation of the SEIV in different energy intervals can be 
up to 1.37x10? times for producing Cf-252. Therefore, the efficiency of 
transplutonium isotopes can be improved by modifying the neutron spectrum around 
the target [42-43]. 

As shown in Table 4, the ESTV varies greatly for different irradiation schemes. 
For example, the maximum deviation of the ESTV can be up to 21.3 times for 
producing Cf-252. Moreover, the bare target always has the lowest ESTV regardless 
of the isotope produced. So the production of transplutonium isotopes requires placing 
a moderator around the target. Layout #6 of the moderator has the highest ESTV for 
producing Pu-238, Cm-242, Cm-244, and Layout £5 has the highest ESTV for 
producing Cf-252. Since the SEIV and ESTV can be used to analyze the production 
efficiency of each energy interval in detail, the SEIV and ESTV are helpful in 
providing direction for optimizing the irradiation scheme. 


Table 3. SEIV for producing Pu-238, Cm-242, Cm244 and Cf-252 


Energy / MeV Pu-238 Cm-242 Cm-244 Cf-252 

[1.00E-11, 1.00E- 2.42E+00 
07] 4.55E+04 3.06E+03 4.96E+00 

[1.00E-07, 4.14E- 2.53E+00 
07] 2.78E+04 2.10E+03 1.10E+01 

[4.14E-07, 8.76E- 1.75E+01 
07] 6.99E+04 1.88E+03 9.02E-01 

[8.76E-07, 1.86E- 1.16E+02 
06] 1.93E+04 4.87E+02 7.34E-02 

[1.86E-06, 5.04E- 3.89E+03 3.17E+02 1.89E-02 1.48E+00 
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06] 
[5.04E-06, 1.07E- 
05] 
[1.07E-05, 3.73E- 
05] 
[3.73E-05, 1.01E- 
04] 
[1.01E-04, 2.14E- 
04] 
[2.14E-04, 4.54E- 
04] 
[4.54E-04, 1.58E- 
03] 
[1.58E-03, 3.35E- 
03] 
[3.35E-03, 7.10E- 
03] 
[7.10E-03, 1.50E- 
02] 
[1.50E-02, 2.19E- 
02] 
[2.19E-02, 2.42E- 
02] 
[2.42E-02, 2.61E- 
02] 
[2.61E-02, 3.18E- 
02] 
[3.18E-02, 4.09E- 
02] 
[4.09E-02, 6.74E- 
02] 
[6.74E-02, 1.11E- 
01] 

[1.1 LE-01, 1.83E- 
01] 
[1.83E-01, 2.97E- 
01] 
[2.97E-01, 3.69E- 
01] 
[3.69E-01, 4.98E- 
01] 
[4.98E-01, 6.08E- 
01] 
[6.08E-01, 7.43E- 
01] 
[7.43E-01, 8.21E- 
01] 
[8.21E-01, 
1.00E+00] 
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7.91E-03 


4.01E-03 


1.17E-01 


2.15E-01 


1.52E-01 


1.70E-01 


2.62E-01 


1.16E-01 


1.19E-01 


8.62E-02 


3.31E-02 


1.62E-02 


1.36E-02 


].15E-02 


9.06E-03 


6.17E-03 


3.88E-03 


2.01E-03 


1.48E-03 


1.19E-03 


8.20E-04 


4.86E-04 


2.75E-04 


1.48E-04 


6.30E-05 


2.94E-05 


9.41E+00 


4.46E+00 


1.91E+00 


1.90E+00 


1.90E+00 


1.73E+00 


7.63E-01 


5.65E-01 


4.08E-01 


3.14E-01 


2.83E-01 


2.67E-01 


2.46E-01 


2.12E-01 


1.49E-01 


7.85E-02 


4.19E-02 


2.19E-02 


1.34E-02 


7.76E-03 


3.71E-03 


1.85E-03 


1.12E-03 


7.19E-04 


[1.00E+00, 
1.35E+00] 
[1.35E+00, 
1.65E+00] 
[1.65E+00, 
1.92E+00] 
[1.92E+00, 
2.23E+00] 
[2.23E+00, 
2.35E+00] 
[2.35E+00, 
2.37E+00] 
[2.37E+00, 
2.47E+00] 
[2.47E+00, 
2.73E+00] 
[2.73E+00, 
3.01E+00] 
[3.01E+00, 
3.68E+00] 
[3.68E+00, 
4.97E+00] 
[4.97E+00, 
6.07E--00] 
[6.07E--00, 
7.41E+00] 
[7.41E--00, 
8.61E+00] 
[8.61E+00, 
1.00E+01] 
[1.00E+01, 
1.22E+01] 
[1.22E+01, 
1.42E+01] 


5.54E-04 


2.45E-04 


1.39E-04 


9.06E-05 


6.45E-05 


5.95E-05 


5.53E-05 


4.52E-05 


3.26E-05 


2.3] E-05 


1.18E-05 


5.27E-06 


2.35E-06 


1.17E-06 


7.86E-07 


5.51E-07 


0 


1.45E-03 


6.23E-04 


3.14E-04 


1.39E-04 


7.85E-05 


6.54E-05 


5.63E-05 


3.94E-05 


2.00E-05 


7.95E-06 


1.17E-06 


1.05E-07 


1.77E-08 


9.91E-09 


1.13E-08 


0.00E+00 


6.32E-09 


1.12E-05 


4.44E-06 


2.50E-06 


1.43E-06 


9.42E-07 


8.18E-07 


7.20E-07 


5.13E-07 


3.04E-07 


1.85E-07 


1.10E-07 


3.67E-08 


4.22E-09 


4.80E-09 


2.56E-08 


1.19E-07 


1.32E-07 


4.20E-04 


3.30E-04 


2.82E-04 


1.93E-04 


1.35E-04 


1.20E-04 


1.10E-04 


8.03E-05 


4.89E-05 


2.09E-05 


3.64E-06 


4.23E-07 


1.37E-07 


8.47E-08 


8.52E-08 


8.78E-08 


0 


Table 4. ESTV for producing Pu-238, Cm-242, Cm244 and Cf-252 


Schemes 
Bare target 
Layout#l 
Layout#2 
Layout#3 
Layout#4 


Layout#5 


Pu-238 
4.50E-03 
1.36E-01 
6.78E-01 
1.60E+00 
5.46E+00 
6.26E+00 


Cm-242 
5.64E-03 
1.92E-02 
4.05E-02 
7.1 4E-02 
1.74E-01 


1.99E-01 


Cm-244 
3.92E-06 
2.55E-05 
5.00E-05 
7.59E-05 
1.77E-04 


2.01E-04 


Cf-252 
9.22E-05 
4.41E-04 
1.08E-03 
1.72E-03 
4.39E-03 
4.94E-03 
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Layout#6 1.11E+01 2.94E-01 2.42E-04 4.59E-03 

It is not rigorous to compare the different irradiation schemes only by the SEIV 
and ESTV, so we need to verify this method. To verify the effectiveness of the rapid 
diagnosis method, we perform the burnup calculations with RMC code [44] to 
calculate the yields of Pu-238, Cm-242, Cm-244, and Cf-25, respectively. Since the 
neutron spectrum is changing during the burnup process, we first simulate the one- 
second burnup process to ensure that the neutron spectrum used to calculate the ESTV 
is consistent with the neutron spectrum used to calculate the yields. In the one-second 
burnup process, we can clarify the relationship between the ESTV and the yields with 
a fixed neutron spectrum. The yields in one second for Pu-238, Cm-242, Cm244 and 
Cf-252 are given in Table 5. 


Table 5. Yields in one second for Pu-238, Cm-242, Cm244 and Cf-252 (/g) 


Schemes Pu-238 Cm-242 Cm-244 Cf-252 
Bare target 5.03E-13 ].53E-12 1.68E-45 1.94E-37 
Layout#1 1.07E-12 2.29E-12 4.96E-45 6.90E-37 
Layout#2 1.86E-12 4.20E-12 1.07E-44 1.32E-36 
Layout#3 2.86E-12 6.75E-12 1.92E-44 1.79E-36 
Layout#4 6.78E-12 1.62E-11 2.44E-43 1.38E-35 
Layout#5 7.62E-12 1.80E-11 2.85E-43 1.46E-35 
Layout#6 1.05E-11 2.64E-11 1.60E-43 5.79E-36 


As shown in Table 5, no matter which isotope is produced, the yield in one- 
second increases with the increases of the ESTV. Therefore, the ESTV can be used to 
evaluate the yield with the same neutron spectrum. However, the one-second burnup 
can only describe the physical process at the early stage of irradiation but cannot 
describe the entire irradiation cycle. To ensure whether the ESTV can indicate the 
yield of the whole irradiation cycle, the burnup process of the target irradiation for 90 
days is also simulated. The yields in 90 days for Pu-238, Cm-242, Cm244 and Cf-252 


are given in Table 6. 
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Table 6. Yields in 90 days for Pu-238, Cm-242, Cm244 and Cf-252 (/g) 


Schemes Pu-238 Cm-242 Cm-244 Cf-252 
Bare target 1.65E+00 1.28E+00 1.96E-02 7.95E-03 
Layout#1 3.03E+00 1.74E+00 4.5] E-02 2.07E-02 
Layout#2 4.64E+00 2.68E+00 7.91E-02 3.42E-02 
Layout#3 6.07E+00 3.58E+00 2.98E-01 4.20E-02 
Layout#4 9.44E+00 6.57E+00 1.00E+00 8.23E-02 
Layout#5 1.01E+01 7.08E+00 7.08E+00 9.80E-02 
Layout#6 1.16E+01 1.05E+01 1.05E+01 8.30E-02 


As shown in Table 6, no matter which isotope is produced, the yield in 90 days 
increases with the increase of the ESTV. Therefore, the ESTV can not only indicate 
the production efficiency at the early stage but also indicate the efficiency of the 
whole irradiation cycle. So the rapid diagnosis method based on the SEIV and ESTV 


can be used to evaluate the efficiency of transplutonium isotopes production. 


3.2 Experimental Verification 

Since the high flux reactor shown in Figure 4 has yet to be built, experimental 
studies cannot be carried out on this reactor, so the experimental verification of the 
rapid diagnosis method is carried out on the High Flux Isotope Reactor (HFIR) [45]. 
HFIR has been producing Cf-252 since 1966 and currently accounts for 70% of the 
world’s supply of Cf-252, so HFIR has a large number of measured data in 
transplutonium isotopes production. 


Rated at 100 MW and currently operating at 85 MW, HFIR has the world’s 
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highest steady-state neutron heat flux (2.6x10P cm?:s'). HFIR is a light-water-cooled 
moderated flux trap reactor using highly enriched U-235. The reactor core consists of 
many concentric annular zones, each about 61 cm high (the fuel height is 51 cm). At 
the center is a 12.70 cm cylindrical pore called “neutron flux trap", containing 37 


vertical experimental holes. The neutron flux trap is surrounded by two concentric 


fuel assemblies. The reactor core and neutron flux trap are shown in Figure 10 [46]. 


reactor core neutron flux trap 


Figure 10. The reactor core and neutron flux trap of HFIR 
Because the rapid diagnosis method uses the ESTV to evaluate the production 
efficiency of different irradiation schemes, the Monte Carlo criticality calculation 
should be performed to calculate the ESTV. RMC code is used to model HFIR, and 
the input card of HFIR is given in the attachment. See Figure 11 for the core modeling 


and the neutron flux trap modeling. 


reactor core neutron flux trap 
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Figure 11. The core modeling and the neutron flux trap modeling of RMC code 
To verify the correctness of the modeling, the radial distributions of thermal 
neutron flux and total neutron flux are calculated and compared with reference 
solutions, as shown in Figure 12. It shows that the flux distributions calculated by 


RMC are in good agreement with the reference solutions, which proves the 


correctness of the HRIR modeling. 


a 


6E+13 


5E+13 


=E= RMC-thermal flux 
6 - Ref-thermal flux 
—À-— RMC-total flux 
—W-~ Ref-total flux 


4E+13 


3E+13 


Neutron Flux (cm72/s) per MW Power 


2E+13 


Radial Distance from Reactor Center (cm) 

Figure 12. Comparison of the radial flux distributions 
Then experimental verification is carried out. The target design shown in Figure 
13 was irradiated in HFIR for five days. Four irradiation schemes were formed with 
two target nuclide compositions and two irradiated positions, namely, ??Cm A, 
?5Cm C, Cf A, and Cf C. The yields of heavy nuclei (all nuclides heavier than the 
nuclides in the target) for the four irradiation schemes were measured [9], and the 


corresponding ESTV was calculated. The experimental results are shown in Table 7. 
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Figure 13. Experimental target design 


Table 7. Results of experimental verification 


Irradiation 
schemes 


245Cm A 
245m C 


Cf A 


Cf C 


Nuclide composition Irradiation pA ; ESTV 
i heavy nuclei : 
of the target duration . (simulations) 
(experiments) 
245Cm: 45.9+0.3ng 7.861 ng 1.00E-12 
5 days 
249 S 
249 š 
CE A2 DOE 12.6 ng 323E-11 
CF: 11.340.2ng 5 days 
12.7 ng 3.50E-11 


21Cf: 30.2+0.5ng 


Because the target nuclide composition should be the same when comparing the 


ESTV of different irradiation schemes, we compared 245Cm A with 245Cm C and 


Cf A with Cf C, respectively. As shown in Table 7, the measured yield of heavy 


nuclei of ??Cm A is lower than that of °*⁄ Cm C, and the corresponding ESTV of 


?5Cm A is also lower than that of "Cm C. Meanwhile, the measured yield of heavy 


nuclei of Cf A is lower than that of Cf C, and the corresponding ESTV of the Cf A 


is also lower than that of Cf C. The experimental results satisfy the law in the rapid 
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1 diagnosis method that the higher the ESTV is, the higher the production efficiency is. 
2 Therefore, the experimental verification proves the correctness of the rapid diagnosis 
3 method. 

4 In the next chapter, we will use the rapid diagnosis method to realize the rapid 


5 screening and determination of the optimal irradiation scheme. 
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22 4. SENSITIVITY ANALYSIS AND OPTIMIZATION APPLICATION 

23 4.1 The Logic of Optimization 

24 There are five design parameters to be optimized: 1) the selection of the 
25 moderator, 2) the layout of the moderator, 3) the position of the target, 4) the size of 
26 the target, and 5) the shape of the target. 


27 As shown in Table 3, no matter which isotope is produced, the energy interval 
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with the highest SEIV is always in the thermal range, so we need to arrange the 
moderator to soften the neutron spectrum around the target. We hope to select the 
material with the strongest moderating power as the moderator so that we can achieve 
the best neutron spectrum by deploying the least amount of moderator. The 
moderating power of various materials is the inherent property of the material and is 
not affected by other design parameters, so the moderator material can be determined 
first. We use a large hexagonal prism target with Layout #1 and Position #1 as the 
initial scheme. Nine moderator materials are tested, and the neutron spectrums around 


the target are tallied, as shown in Figure 14. 


10° 
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10° 
104 
10° benzene 
copper 
-6 
% 10 nickel 
197 beryllium 
— —- light water 
10° zirconium hydirde 


stainless steel 
10? graphite 
heavy water 


107 10° #105  10* 10? 10? 10 10? 10! 
Energy / MeV 


Figure 14. Neutron spectrums with different moderators 
As shown in Figure 14, zirconium hydride has the best moderating power, so 
zirconium hydride is selected as the moderator, and the moderator selection will not 
be considered in all subsequent analyses. 
To determine an optimal production scheme, we have to decide the optimal 
design parameter for the layout of the moderator, the position, size and shape of the 


target. We use a large hexagonal prism bare target with Position #1 as the initial 
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scheme and then optimize the four design parameters. The optimization order of these 
four parameters is important. We determine the optimization order by sensitivity 
analysis. The sensitivity coefficient of a design parameter d to production efficiency is 


defined as 
j (10) 


where P represents the production efficiency, d represents a design parameter. 

So the sensitivity coefficient is positively correlated with the ESTV's maximum 
increase that can be obtained by adjusting a design parameter. The larger the 
maximum increase, the more important this design parameter is to the optimization 
process, so we optimize it first. After one of these design parameters is optimized, we 
take the optimized scheme as the starting scheme to optimize the remaining design 
parameters, and so on, until all design parameters are optimized. A logical diagram of 
the optimization process is shown in Figure 15. This paper takes Cf-252 as an 
example to describe the optimization process of the irradiation scheme in detail. The 


optimal schemes for Pu-238, Cm-242 and Cm-244 will be given in Section 3.6. 


Optimization#1 D Size Layout Position Shape 


Optimization&2]) 


_ max(ESTV;) 
- min(ESTV) 


Optimization#3 


N4,1| determined 


l 1] represents the maximum 

increase that the ESTV | 
: can obtain when adjusting : 
Final scheme l this design parameter | 


Optimization£4]) 


determined 


Figure 15. The diagram of optimization logic 


4.2 First Round Sensitivity Analysis and Optimization 
The first round of sensitivity analysis and optimization needs to consider four 


design parameters, i.e., the layout of the moderator, the position, size and shape of the 
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target. The optional values for the four design parameters are shown in Figure 9. New 
irradiation schemes can be obtained by adjusting one of the four design parameters 


based on the initial scheme. The ESTVs for all new schemes are shown in Table 8, 


and the symbol “y” represents the sensitivity coefficient. 

Table 8. The ESTVs for the first round optimization 
Layout ESTV Position ESTV 
bare target 9.22E-05 Position#1 9.22E-05 
Layout#1 4.41E-04 Position#2 1.07E-04 
Layout#2 1.08E-03 Position#3 5.43E-04 
Layout#3 1.72E-03 Position#4 6.97E-04 
Layout#4 4.39E-03 
Layout#5 4.94E-03 
Layout#6 4.59E-03 
N. 53.6 np 7.56 
Shape ESTV Size ESTV 
hexagonal prism 9.22E-05 large 9.22E-05 
sphere 9.78E-05 medium 9.35E-05 
cylinder 9.22E-05 small 9.22E-05 
ysn 1.06 si 1.01 

As shown in Table 8, the layout of the moderator can improve the ESTV the 


most, so the layout of the moderator should be optimized first. Moreover, Layout#5 
has the highest ESTV, so the layout of the moderator is selected as Layout #5. After 
the first round of optimization, the scheme becomes a large hexagonal prism target 
with Position #1 and Layout #5. 
4.3 Second Round Sensitivity Analysis and Optimization 

The second round of sensitivity analysis and optimization needs to consider three 


design parameters, i.e., the position, size and shape of the target. All new irradiation 
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schemes are obtained based on the scheme after the first round of optimization. The 
ESTVs for the second round of optimization are shown in Table 9. 


Table 9. The ESTVs for the second round optimization 


Position ESTV Shape ESTV Size ESTV 


Position#l  4.94E-03 | hexagonal prism ^ 4.94E-03 | Large 4.94E-03 
Position#2  3.45E-03 | sphere 4.03E-03 | Medium 5.42E-03 
Position#3 2.23E-03 | cylinder 4.75E-03 | Small 6.28E-03 
Position#4 1.04E-03 


n» 4.75 Nsh 1.04 si 1:27 

As shown in Table 9, the target's position can improve the ESTV the most, so the 
target's position should be optimized in the second round of optimization. Moreover, 
Positionz1 has the highest ESTV, so the target's position is selected as Position #1. 
After the second round of optimization, the scheme becomes a large hexagonal prism 
target with Position #1 and Layout #5. 
4.4 Third Round Sensitivity Analysis and Optimization 

The third round of sensitivity analysis and optimization needs to consider two 
design parameters, i.e., the size and shape of the target. All new irradiation schemes 
are obtained based on the scheme after the second round of optimization. The ESTVs 
for the third round of optimization are shown in Table 10. 


Table 10. The ESTVs for the third round optimization 


Shape ESTV Size ESTV 
hexagonal prism 4.94E-03 Large 4.94E-03 
sphere 4.03E-03 Medium 5.42E-03 
cylinder 4.75E-03 Small 6.28E-03 
Nsh 1.04 si 1.27 


As shown in Table 10, the target's shape and size have similar effects on 
improving the ESTV, with the target’s size slightly larger, so the target’s size should 


be optimized in the third round optimization. Moreover, the small target has the 
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highest ESTV, so the small target is selected. After the third round of optimization, the 
scheme becomes a small hexagonal prism target with Position #1 and Layout #5. 


4.5 Fourth Round Sensitivity Analysis and Optimization 


The fourth round of sensitivity analysis and optimization needs to consider one 
design parameter, i.e., the shape of the target. All new irradiation schemes are 
obtained based on the scheme after the third round of optimization. The ESTVs for 
the fourth round of optimization are shown in Table 11. 

Table 11. The ESTVs for the fourth round optimization 
Shape ESTV 

hexagonal prism 6.28E-03 

sphere 5.19E-03 

cylinder 6.10E-03 

Nsn 1.21 

As shown in Table 11, the hexagonal prism has the highest ESTV, so a hexagonal 
prism target is selected. After the fourth round of optimization, the scheme becomes a 
small hexagonal prism target with Position #1 and Layout #5. 

4.6 Final Irradiation Scheme 

The optimized scheme for producing Cf-252 is a small hexagonal prism target 
with Position #1 and Layout #5. Similarly, we optimized the schemes for producing 
Pu-238, Cm-242 and Cm-244. The final optimized irradiation schemes for producing 
Pu-238, Cm-242, Cm-244 and Cf-252 are given in Table 12. The comparison of 
production efficiency before and after optimization is shown in Figure 16. 

Table 12. The final optimized irradiation schemes 
Design parameters Pu-238 Cm-242 Cm-244 Cf-252 
zirconium zirconium zirconium zirconium 
Selection of moderator hydride hydride hydride hydride 


Layout of moderator Layout #6 Layout #6 Layout #6 Layout #5 


Position of target Position #1 Position #1 Position #1 ^ Position #1 
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Shape of target hexagonal hexagonal cylinder hexagonal 


Size of target small small small small 


initial scheme for Cm-242 
optimized scheme for Cm-24 


initial scheme for Pu-238 
optimized scheme for Pu-238 
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Figure 16. The comparison of the efficiency before and after optimization 


After four rounds of sensitivity analysis and optimization, we can quickly obtain 
the optimal scheme for transplutonium isotopes production. As shown in Figure 16, 
the optimal scheme can effectively improve production efficiency compared with the 
initial scheme. The efficiency of Pu-238 is increased by 7.41 times, 11.98 times for 
Cm-242, 65.20 times for Cm-244 and 15.08 times for Cf-252. Therefore, the 
optimization strategy based on the rapid diagnosis method is helpful for 


transplutonium isotopes production. 
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5. CONCLUSION 


The production of transplutonium isotopes with reactor irradiation is inefficient, 
and optimization studies are needed. The physics process of transplutonium isotopes 
production is complicated to analyze quantitatively. After grasping the physics nature 
and simplifying the complexity of the production process, we proposed the concept of 
“Single Energy Interval Value (SEIV)" and “Energy Spectrum Total Value (ESTV)” to 
realize the quantitative analysis of the production efficiency of transplutonium 
isotopes. A rapid diagnosis method for evaluating the radiation scheme based on the 
SEIV and ESTV is established. The optimization design for producing Cf-252, Cm- 
244, Cm-242 and Pu-238 is carried out based on the rapid diagnosis method. 
Compared with the initial scheme, the optimal scheme improves the efficiency of Pu- 
238 by 7.41 times, 11.98 times for Cm-242, 65.20 times for Cm-244 and 15.08 times 
for Cf-252. 

This paper provides solutions to two key questions in the production of 
transplutonium isotopes: 1) how to find the optimal neutron spectrum for producing 
transplutonium isotopes, and 2) how to achieve the optimal neutron spectrum. For the 


first question, we use the SEIV and ESTV to achieve a detailed analysis of the 
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neutron spectrum around the target. For the second problem, we optimized the design 


parameters around the target to achieve the optimal neutron spectrum, effectively 


improving the computational efficiency. 


The work in this paper realizes the refined analysis of transplutonium isotopes 


production and provides a theoretical basis for improving production efficiency. We 


will research energy spectrum conversion technology to achieve the optimal energy 


spectrum more finely in the future. 
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